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ABSTRACT 

We present the J and H —band source catalog covering the AKARI North Ecliptic Pole field. Filling 
the gap between the optical data from other follow-up observations and mid-infrared (MIR) data 
from AKARI , our near-infrared (NIR) data provides contiguous wavelength coverage from optical to 
MIR. For the J and FI—band imaging, we used the FLoridA Multi-object Imaging Near-ir Grism 
Observational Spectrometer (FLAMINGOS) on the Kitt Peak National Observatory 2.1m telescope 
covering a 5.1 deg 2 area down to a 5<r depth of ~21.6 mag and ~21.3 mag (AB) for J and H —band 
with an astrometric accuracy of 0.14" and 0.17" for ler in R.A. and Deck directions, respectively. We 
detected 208,020 sources for J—band and 203,832 sources for IF—band. This NIR data is being used 
for studies including analysis of the physical properties of infrared sources such as stellar mass and 
photometric redshifts, and will be a valuable dataset for various future missions. 
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1. INTRODUCTION 

AKA RI is an infrared (IR) space telescope launched 
in 2006 (lMurakami et al.l 1200 7) which has conducted an 
all-sky survey in far-IR (FIR) as well as imaging obser¬ 
vations of regions of special interest from near-IR (NIR) 
to FIR. Since AKARI was placed on a Sun-synchronous 
orbit, it conducted deep observations only toward the 
north and the south ecliptic poles. Observing these areas 
in the contiguous nine passbands between 2.4 and 26/iin 
covering NIR to mid-IR (MIR), AKARI provided valu¬ 
able NIR to MIR imaging survey data where the Spitzer 
Space Telescope has limitations, especially at the 11, 15, 
and 18/rm bands. 

Taking the advantage of the Sun-synchronous orbit, 
AKARI carried out deep surveys with nine filters of the 
InfraRed Camera (IRC) that cover from 2 to 26 /im 
contiguously t oward the region around the North Eclip¬ 
tic Pole fNEP: IMatsuhara et al.ll2006t iWada et al.ll2007t 
iLee et all 12007112009D . The NEP survey is composed of 
two main programs, N EP-Wide and NEP-Deep surveys . 
The NEP-Wide survey (ILee et, al.ll2009l iKim et all 1 20 121) 
covers the entire 5.4 deg 2 to the depths of 20.9 (AB mag) 
to 17.8 (AB mag) from 2.4 to 24 p m (5 a point-source 
sensitivity), whereas the NEP-Deep survey (IWada et al.1 
120081: Ffakagi et, all [Mil covers a 0.67 deg 2 area with 
exposures 3-10 times longer than the NEP-Wide sur- 
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vey. With its continuous MIR wavelength range and 
wide field of view, AKARI NEP survey provides a unique 
dataset for the study of galaxy evolution, and it has been 
used to unveil the link between the star formation and 
the active galactic nucleus activity of dis tant galaxies 
(jKarouzos et al.1 12014! lHanami et al.l 1 201 2D, the cosm ic 
star formation history beyond z > 1 dGoto et al] [2010'). 
the properties of polycyclic aroma t ic hydrocarbons him i- 
nous galaxies (iTakagi et al.l 1201(1 lOhvama et al.1 1200711 . 
the environmental dependence of star formation activ¬ 
ities of galaxies in a superc l uster and MIR-excess of 
early-type galaxies dKo et al.ll2012D . an d the nature of 
sources in 11 and 15 pm . limited samples ([Lee et al.ll2007t 
IMatsuhara et al.l 120071 : iPearson et all 1201011 Derivation 
of the accurate local luminosity functions in MIR is un¬ 
der progress (Kim et al., in preparation) , by utilizing the 
subsequent spectroscopic survey data dShim et al .1120131 . 

For the AKARI NEP field, a considerable amount 
of ground-based ancillary datasets have been obtained. 
Figure [T] shows the filter response curves of the multi¬ 
wavelength datasets used in the AKARI NEP field sur¬ 
vey. Follow up imaging data observed by the Mega- 
Cam of the Canada-France-Hawaii Telescope (CFHT) 
with and z’ filters covers 2 deg 2 of th e cen¬ 

tral region of the NEP field dHwang et, all 120071 thick 
dotted blue lines in Figure [lj. Also for the remaining 
area, I Jeon et all ( 2010 ) provided optical data of B , R , 
and I—band covering 4.9 deg 2 (thick gray lines in Figure 
[I]) using the Seoul Nati onal JJniversity 4k x 4k Cam¬ 
era (SNUCAM; Urn et, all 120101) on the 1.5m Ritchey- 
Chretian, AZT-22 telescope at the Maidanak observa¬ 
tory. Spectroscopic data of 1796 sources were obtained 
using Hectospec of MMT and Hydra of the WIYN tele¬ 
scope (|Shim et al.1 12013D . The radio data were taken 
with the Westerbork Radio Synthesi s Telescope (WSRT ) 
over a 1.7 deg 2 area around NEP (IWliite et al.l l2010h . 
and ultraviolet (UV) imaging data were obtained with 
the Galaxy Evolution Explorer (GALEX) at Far-UV and 
Near-UV (Malkan et al, in preparation). More re¬ 
cently, the entire A7C4AINEP field has been covered by 
the Spectral and Photometric Imaging Receiver (SPIRE; 
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Figure 1. Filter response curves of multi-wavelength data covering the NEP field. The thick dotted blue lines represent the 
and z' band filters from MegaCam of CFHT. The thick gray lines are for the B, R, and /—band filters from Maidanak observatory. The 
thin solid red lines and thin dashed red lines are for the AKARI IR. band filter coverage (TV2, N 3, NA , S 7, S9W, 511, L15, Ll&W, and L24). 
The thick black lines are for our J and H —band filters, which make the wavelength coverage between the optical data and the AKARI 
data continuous. Table [l] lists the information of multi-wavelength data for the filters, shown here. 


iGriffin et al.ll2010l) and a part of the NEP area is covered 
by the Photoconductor Array Camera and Spectrometer 
(TAG'S: IPoglitsch et alJl20 1(l) o f the Herschel Space Ob¬ 
servatory dSerieant et, al.ll2012li . In addition, the NEP- 
Deep area has been observed for optical/NIR. im ages with 
the Subaru telescope and the CFHT (lOi et al.ll2014f) and 
in optical spectroscopy with the DEep Imaging Multi- 
Object Spectrograph (DEIMOS) on the Keck telescope 
(Takagi et al., in preparation). We summarized the infor¬ 
mation of the multi-wavelength surveys for the optical, 
NIR and MIR bands in Table [1] for NEP-Wide and Table 
[2] for NEP-Deep. _ _ 

The AKARI data (iKirn et all 120121 thin solid red 
lines and thin dashed red lines in Figure |TJ) covers a 
wavelength range from 2.4 to 24/xm with nine filters 
{N2,N3,N4,S7,S9W,SU,L15,L18W, and L24, where 
the number indicates the central wavelength of the filter 
in pan), while the ground-based imaging data offering 
the coverage at A < 0.9 pm. However there is a gap 
in the wavelength coverage between the optical /—band 
(0.85pm) and the shortest AKARI band (2.4pm). Fill¬ 
ing this gap in the spectral energy distribution (SED) 
will significantly enhance the value of the NEP dataset 
in many ways. First, by adding the NIR data, char¬ 
acterization of the SEDs with a contiguous wavelength 
range from NIR to MIR will increase the accuracy in 
estimating stellar masses of low redshift galaxies, be¬ 
cause the shorter (or longer) wavelength data may be 
contaminated by star formation and dust content. Sec¬ 
ond, the NIR coverage can help improve the photomet¬ 
ric redshift estimation by sampling crucial wavelengths 
at which some prominent features, such as the 1.6pnr 
bump or the 4000A break, may reside for low redshift or 
high redshift galaxies, respectively. 

The NEP field has served as a major survey field 
for extragalactic astronomy for AKARI , and Herschel 
missions, but it is also poised to do so in future mis¬ 
sions such as Euclid (jLaureijs et al.ll20TTI) , Space Infrared 
Telescope for Cosm ology and Astrophysics ( SPICA ; 
iSerieant et al.1 12012), Cosmic Infrare d Background Ex¬ 
periment ( CIBER ; Bock et all I2013T) . extended ROent¬ 


gen S urvey with an Imaging Telescope Array (eROSITA; 
iMerloni et al.ll2 012fl. and LOw-Frequency ARray (LO- 
FAR: Ivan Haarlem et all 1 201, ‘11 1. In order to fill the gap 
in the wavelength coverage and to enhance the scien¬ 
tific value of the AKARI NEP survey, we conducted an 
imaging observation of the entire NEP survey area in J, 
and H- band (thick black lines in Figure [l]). Here, we 
present the J and H— band data of the AKARI NEP 
field. This paper is organized as follows. In Section O 
we describe the details of the observations. In Section 
03 data reduction processes are explained. In Section [2 
we describe the source detection and photometry along 
with the properties of the data and provide the J and 
//—band merged catalog. In Section 0 we provide the 
result from our data such as source number counts and 
color-magnitude diagrams. We summarize our results 
in the final section. Throughout this paper, we use a 
cosmology with Hm = 0.3, Ha = 0.7, and H 0 = 70 km 
s _1 Mpc^ 1 . We use the AB magnitude system. 

2. OBSERVATIONS 

The AKARI NEP held is centered at a = 
18 h 00 m 00 s , 5 = +66°33'38" and our NIR data cover the 
entire AKARI ’s NEP held area. Figure [2] (left) shows 
tiling components of our NIR data (small empty green 
squares; each of them has a 21' x 21' held of view.) as 
well as the held coverage of the AKARI NEP surveys 
(solid red circles) and its multi-wavelength follow-up ob¬ 
servations. The two dashed blue squares show the cover¬ 
age of the CFHT observation (u*g'r’i'z'— bands) and the 
gray area is for the Maidanak observation (BRI— bands). 
The four filled orange squ ares have J and A%—band 
imaging data introduced in llmai et all (120071 1 covering 
only the NEP-Deep area. The right panel of Figure [2] 
shows the same held coverage as the left one, but shows 
the held name for each subheld of our NIR data. From 
FLM01 to FLM66, 65 tiles constitute the NEP NIR area. 
Three subhelds, FLM64 to FLM66, have only //—band 
data from our o bserv ation, because this area overlaps 
with llmai et all (|2007f) ’s J and K s —band data. The tile 
names are used for the target ID (Section 14.51) . The 
names of each subheld and their central positions are 
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Table 1 

Information of the multi-wavelength survey for the optical, NIR and MIR bands in NEP-Wide survey. 


Instrument 

Area 

Pixel Size 

Filter 

Effective Wavelength 

Sensitivity 

Reference 


(deg 2 ) 

(" Pix" 1 ) 


H 

(5cr, AB mag) 


MegaCam (CFHT) 

2 

0.19 

u* 

0.37 

25.8 

Hwang et al. (2007) 




a' 

0.49 

25.9 




r' 

0.63 

25.4 





i' 

0.77 

24.6 





z' 

0.89 

23.7 


SNUCAM (Maidanak 1.5m) 

4.9 

0.27 

B 

0.44 

23.4 

Jeon et al. (2010) 




R 

0.64 

23.1 





I 

0.79 

22.3 


FLAMINGOS (KPNO 4m) 

5.1 

0.30 

J 

1.25 

21.6 

This work 



H 

1.64 

21.3 


IRC (AKARI) 

5.4 

1.46x1.46 

N2 

2.4 

20.9 

Kim et al. (2012) 


N3 

3.2 

21.1 

N4 

4.1 

21.1 

2.34X2.34 57 

7.0 

1975 

S9W 

9.0 

19.3 

S 11 

11.0 

19.0 

2.51x2.39 L15 

I3D 

18.6 

L18W 

18.0 

18.7 

L24 

24.0 

17.8 


Table 2 

Same as Table [l] but in NEP-Deep survey. 

Instrument 

Area 

(deg 2 ) 

Pixel Size 
(" Pix -1 ) 

Filter Effective Wavelength 

Om) 

Sensitivity 
(5cr, AB mag) 

Reference 

MegaCam (CFHT) 

1 

0.19 


0.37 

24.6 

Takagi et al. (2012) 




a' 

0.49 

26.5 

Oi et J. (2014) 




r' 

0.63 

25.7 





i' 

0.77 

24.9 





z' 

0.89 

23.9 


SuprimeCam (Subaru) 

0.255 

0.2 

B 

0.43 

28.4 

Takagi et al. (2012) 




V 

0.54 

28.0 





R 

0.65 

27.4 





i' 

0.80 

27.0 





z' 

0.91 

26.2 


WIRCam (CFHT) 

0.6 

0.3 

Y 

1.02 

23.2 

Oi et al. (2014) 




J 

1.25 

22.8 





I<s 

2.15 

22.5 


FLAMINGOS (KPNO 4m) 

0.2 

0.61 

J 

1.25 

21.9 

Imai et al. (2007) 




K s 

2.15 

21.0 


IRC (AKARI) 

0.67 

1.46x1.46 

N2 

2.4 

21.5 

Takagi et al. (2012) 




N3 

3.2 

21.7 





N4 

4.1 

22.1 




2.34x2.34 

S 7 

7.0 

1577 





S9W 

9.0 

19.5 





S 11 

11.0 

19.3 




2.51x2.39 

L15 

reii 

1577 





L1SW 

18.0 

18.7 





L24 

24.0 

17.8 



listed in Table [3] 

The NIR imaging observation was done with J and 
H filters from 2008 June 13 to 23 at Kitt Peak Na¬ 
tional Observatory (KPNO) using the FLoridA Multi¬ 
object Imaging Nea r-ir Grisrn O bserva tional Spectrome¬ 
ter (FLAMINGOS: lElston et alJl2006 ). which is an NIR 
camera attached to the Cassegrain focal plane of the 2.1m 
telescope. It uses a Hawaii II 2048x2048 HgCdTe array, 
divided into four quadrants with 8 amplifiers each. It 
provides a 21' x 21' field of view, with a pixel scale of 
0.606" pix” 1 . 

The observation was carried out under mostly clear 
sky conditions. The observations followed a standard se¬ 
quence which consists of a 5 x 5 dither pattern with 20" 
offsets to make the sky frame for the sky subtraction. 


Each subfield overlaps with neighboring subfields for re¬ 
gions of width 4' to avoid the low efficiency parts on the 
edges of each chip. The total exposure time for each 
subfield is ~36 minutes for J—band and ~21 minutes for 
H— band, and the data reach 5er depths of ^21.6 mag and 
~21.3 mag, respectively. We used 1.5 minutes as an ex¬ 
posure time per frame for J—band throughout the obser¬ 
vations. However, for the H —band observations, we em¬ 
ployed exposure times between 0.5 and 1.0 minutes per 
frame to adjust the background level. The seeing values 
for each subfield were determined as the median values of 
full width at half maximum (FWHM) of non-saturated 
point sources with magnitudes between 15 and 19. The 
seeing values of the stacked J and H —band images are 
determined to be ~1.8" and ~1.7", respectively. To de- 
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Figure 2. Field coverage of the AKARI NEP multi-wavelength data. The center is the NEP whose coordinates are a = 18™00 m 00 s ,5 = 
+66°33 / 38 // . Our J and H —band survey (5.1 deg 2 ) is indicated with small empty green squares on both panels (65 subfields). The names 
of the subfields are marked with numbers on the right panel. The names and central positions for each subfield are listed in Table [3] The 
solid red circles and the two dashed blue squares indicate AKARI (5.4 deg 2 ) and CFHT MegaCam (2 deg 2 ) coverage, respectively. The 
gray area represents the BRI imaging survey (4.9 deg 2 ) from Maidanak observatory. The thick red cross indicates the center of NEP-Deep 
field (0.67 deg 2 ). The J and K s band coverage from [Tmai et al.l 1120070 (750 arcmin 2 ) are marked with four filled orange squares centered 
at the NEP-Deep field. 


Table 3 

Names of each subfield in Figure [I] and their central positions. 


Field 

R.A. 

Decl. 

Field 

R.A. 

Decl. 

FLM01 

18:04:18.09 

+67:32:00.71 

FLM66 

17:56:29.34 

+66:28:38.50 

FLM02 

18:01:30.29 

+67:36:38.52 

FLM34 

17:53:58.87 

+66:28:38.50 

FLM03 

17:58:39.75 

+67:36:38.52 

FLM35 

17:51:28.40 

+66:28:38.50 

FLM04 

17:55:49.22 

+67:36:38.52 

FLM36 

17:48:57.93 

+66:28:38.50 

FLM05 

18:07:08.92 

+67:15:36.93 

FLM37 

18:11:22.13 

+66:11:38.51 

FLM06 

18:04:18.63 

+67:16:45.62 

FLM38 

18:08:51.66 

+66:11:38.51 

FLM07 

18:01:30.29 

+67:19:38.53 

FLM39 

18:06:21.19 

+66:11:38.51 

FLM08 

17:58:39.75 

+67:19:38.53 

FLM40 

18:03:48.46 

+66:10:45.92 

FLM09 

17:55:49.22 

+67:19:38.53 

FLM41 

18:01:18.11 

+66:10:23.74 

FLM10 

17:52:58.68 

+67:19:38.53 

FLM42 

17:58:49.78 

+66:11:38.51 

FLM11 

18:10:39.67 

+66:59:10.18 

FLM43 

17:56:19.31 

+66:11:38.51 

FLM12 

18:07:51.48 

+67:02:38.51 

FLM44 

17:53:48.84 

+66:11:38.51 

FLM13 

18:05:00.95 

+67:02:38.51 

FLM45 

17:51:18.37 

+66:11:38.51 

FLM14 

18:02:10.41 

+67:02:38.51 

FLM46 

17:48:47.90 

+66:11:38.51 

FLM15 

17:59:19.88 

+67:02:38.51 

FLM47 

18:10:11.92 

+65:54:38.52 

FLM16 

17:56:29.34 

+67:02:38.51 

FLM48 

18:07:21.38 

+65:54:38.52 

FLM17 

17:53:38.81 

+67:02:38.51 

FLM49 

18:04:30.85 

+65:54:38.52 

FLM18 

17:50:48.27 

+67:02:38.51 

FLM50 

18:01:40.31 

+65:54:38.52 

FLM19 

18:11:49.59 

+66:43:31.48 

FLM51 

17:58:49.78 

+65:54:38.52 

FLM20 

18:09:21.76 

+66:45:38.52 

FLM52 

17:55:59.25 

+65:54:38.52 

FLM21 

18:06:51.28 

+66:45:38.52 

FLM53 

17:53:08.72 

+65:54:38.52 

FLM22 

18:04:20.82 

+66:45:38.52 

FLM54 

17:50:18.18 

+65:54:38.52 

FLM23 

18:01:50.35 

+66:45:38.52 

FLM55 

18:07:11.35 

+65:37:38.53 

FLM24 

17:59:19.88 

+66:45:38.52 

FLM56 

18:04:20.82 

+65:37:38.53 

FLM64 

17:56:49.40 

+66:45:38.52 

FLM57 

18:01:30.29 

+65:37:38.53 

FLM65 

17:54:18.93 

+66:45:38.52 

FLM58 

17:58:39.75 

+65:37:38.53 

FLM26 

17:51:48.46 

+66:45:38.52 

FLM59 

17:55:49.22 

+65:37:38.53 

FLM27 

17:49:18.00 

+66:45:38.52 

FLM60 

17:52:58.68 

+65:37:38.53 

FLM28 

18:11:32.17 

+66:28:38.50 

FLM61 

18:03:00.56 

+65:20:38.51 

FLM29 

18:09:01.69 

+66:28:38.50 

FLM62 

18:00:10.03 

+65:20:38.51 

FLM30 

18:06:31.22 

+66:28:38.50 

FLM63 

17:57:19.50 

+65:20:38.51 

FLM31 

18:03:58.38 

+66:27:09.87 





Decl. [degree] 
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termine the limiting magnitude, we used the aperture 
size corresponding to a diameter of 3 times the FWHM 
following the definition introduced in Section I4.fl Ta¬ 
bles [5] and El show the observation summary including 
the total on-source exposure time, the exposure time per 
frame, the median seeing size, 5a detection limit, and 
50% completeness magnitude (see Section 14.61) for each 
filter. We plot the values of seeing size and 5cr detection 
limits for all subfields for both filters in Figure [3] The 
figure can be used to check the subheld variance of the 
values. Also to determine the quality of the images, we 
show plots of the magnitudes for the 5a detection limits 
(crosses), 50% com plete ness (see Section IXIH) . 99% reli¬ 
ability (see Sect ion 1X71) and photometric zeropoints Z p 
(see Section IX3l) for each subheld in Figure [4] The total 
exposure times for each subheld are plotted on the right 
axis. 

3. DATA REDUCTION 

3.1. Preprocessing 

For the data reduction, we hrst carried out i) linear¬ 
ity correction, ii) dark subtraction and hat helding, iii) 
distortion correction, iv) sky subtraction, v) bad pixel 
and cosmic ray masking, and vi) image alignment and 
combining. All these steps are performed using IRAF 6 . 

To correct the nonlinearity of the NIR detector, we 
used the IRAF irlincor task. For FLAMINGOS, the 
linearity correction can be expressed as a third-order 
polynomial 7 : 

y = x x (a + b x (x/32767) + c x (a;/32767) 2 ), 

where x is the raw signal, y is the linearity corrected 
signal in ADU, and the coefficients are a = 1.00425, b = 
-0.03323, c= 0.04489. 

We used a standard method for dark subtraction and 
hat helding with the IRAF noao. imred. ccdred pack¬ 
age. Master dark frames were created by combining 
dark frame image for each night and exposure time. We 
subtracted the master dark image from each science im¬ 
age. These dark subtracted science images, typically 250 
frames with a background count of 3000 DNs, were com¬ 
bined in median after scaling them to a reference image 
by the mode pixel value to construct master sky hat im¬ 
ages for each night. We compared this master sky hat 
image with hat frames taken from dome dimming lights 
and found that the master sky hat can remove ring struc¬ 
tures in the center of science images more efficiently. The 
hat helding was done by dividing dark subtracted im¬ 
ages by the normalized master sky hat image for each 
night. We removed satellite tracks on the images using 
the IRAF satzap task. 

The distortion correction 8 was performed using a dis¬ 
tortion map constructed from pinhole data. We used 
the IRAF geotran task for the correction. Among the 
various distortion maps, we used the most resent one, 
KP2m04.map, obtained in July 2003. The distortion 
map doubles the number of pixels in each axis and sets 

6 IRAF is distributed by the National Optical Astronomy Ob¬ 
servatory, which is operated by the Association of Universities for 
Research in Astronomy (AURA) under cooperative agreement with 
the National Science Foundation. 

' http://www.noao.edu/kpno/manuals/flmn/flmn.user.html 

8 http://flamingos.astro.ufl.edu/DistortionFiles/index.html 


the pixel scale to 0.303" pix -1 . We compared the astro- 
metric accuracy before and after the distortion correction 
(see Section 13.21 for more details about the astrometry). 
We used the FLM14 J—band field as a test held. The 
astrometric accuracy was improved by up to ~5% after 
the distortion correction. 

To subtract the sky, mask bad pixels and cosmic rays, 
and align and combine the sky subtracted individual im¬ 
ages, we used the Experi men tal Deep Infrared Mosaic- 
ing Software (XDIMSUM: [Stanford et al.lll995D package of 
IRAF. Six adjacent frames were used for making the sky 
image to subtract the sky. The mask images were cre¬ 
ated using the master dark images and master sky hat 
images for each night, excluding bad pixels which have 
unusual values compared with the normal distribution 
of the pixel values. Cosmic rays are removed with a 3cr 
clipping method while combining the images. 

3.2. Astrometry 

We employed SCAMP (lBertinll2006lj which hnds astro- 
metric solution by cross matching between positions of 
input catalog created from a science image without as¬ 
trometry calibration and a reference catalog containing a 
set of accurate R.A. and Deck values. The Two Micron 
All Sky Survey Point Source Catalog (2MASS PSC) was 
used as the reference catalog. The astrometric solution 
was derived for each subheld and each hlter separately. 
We tried to hnd the best parameter value for the degree 
of the polynomial for distortion correction between 3, 4, 
5, and 6 and found that a 5 tft, -order polynomial gave the 
best astrometric accuracy among other options. There¬ 
fore we set the degree of the polynomial to 5. 

To check the astrometric accuracy, 2MASS point 
sources with la position uncertainty less than 0.2" are 
compared with our catalog. Figure[5](a) shows the distri¬ 
butions of the positional differences between our sources 
and the corresponding 2MASS sources in R.A. and Deck 
directions, and (b) and (c) show the histogram of the 
position differences in the two directions, respectively. 
We ht the histograms with a Gaussian function (dotted 
lines) and take the ler value of the best fitted Gaussian 
function as the ler astrometric errors in R.A. and Deck 
direction. Here, the ler errors are around 0.16" for both 
directions. Combining the rms scatter in both directions 
in quadrature, we get the absolute positional accuracy 
of 0.17". Considering that the astrometry error of the 
2MASS sources is about 0.2", our astrometric accuracy 
is quite good, and most of the uncertainty in the astrom¬ 
etry seems to originate from the 2MASS catalog. 

3.3. Photometric Calibration 

For the photometric calibration, we also used stellar 
objects in the 2MASS PSC. We restricted the photomet¬ 
ric reference stars to have the color —0.9 < J-imass — 
K 2 MASS < 0.1 (0 < J' 2 M ASS ~ K 2 mASS < 1 for the Vega 
system; see Section lX5l) . photometric error less than 0.05 
mag, and magnitude fainter than 11.9 mag (11 mag for 
the Vega system). For each field, we derived the differ¬ 
ence between the 2MASS and instrumental magnitudes 
of each source in the field and computed the average 
value of the difference to define the zeropoint (Z p ). The 
zeropoint error was derived from the standard deviation 
of the magnitude offset. The average zeropoint error is 
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Figure 3. Seeing and 5<x detection limits for each subfield. (Triangles for J and squares for H— band) 
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Figure 4. Several m easur es of the image quality and depth of each su bfiel d, such as the 5 a detection limit (crosses), the 50% co mplet eness 
(asterisks; see Section 14.61) , the 99% reliability (diamonds; see Section I4.7D . the photometric zeropoint Z p (triangles; see Section |3. 3 L and 
the total exposure time (squares). 
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Figure 5. (a): Distribution of the differences in position between our sources and the corresponding 2MASS sources in R.A. and Decl. 
directions, (b) and (c): Histograms of the distributions in R.A. and Decl. directions, respectively. The dotted lines indicate the best fitted 
Gaussian function for the histograms. From the Gaussian function, the lcr astrometry errors are 0.14" for R.A. and 0.17" for Decl.. 
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Figure 6. Photometric difference between our sources and the corresponding 2MASS sources as a function of magnitude to check the 
photometric accuracy. Although the systematic offset is very small, the bright sources show average offsets of '--'0.05 for both filters. The 
scatter of the magnitude difference is larger for fainter sources, which may be caused by the photometric uncertainties for the sources. The 
dispersion of the magnitude difference are '--'0.05 mag. These values are corresponding with the zeropoint uncertainty. 
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^0.05 mag for both J and H— band. The Z p for each sub¬ 
field is listed in Table ED for the purpose of checking the 
depth of each field. The Z v is used for converting the DN 
per exposure time to magnitudes, i.e. M = Z p — 2.5 log 
(. DN/sec ). Figure 0] shows the distribution of several 
depth indicators for each field. The depths are checked 
using t he 5o ~ detection limit, the 50% complet enes s (see 
Section l4~6ll . the 99% reliability (see Section l4~7l) . and 
the Z p (triangles), as well as the total exposure time. 

To check the photometric accuracy, we compare the 
magnitudes of 2MASS sources with our sources. Figure 
ED shows the magnitude differences between our sources 
and the corresponding 2MASS sources as a function of 
magnitude for J (left panel) and H (right panel) band. 
We used the 2MASS point sources of same color and pho¬ 
tometric error properties that we used for the Z p calcula¬ 
tion. The magnitude offsets show very small systematic 
differences (0.002 for J—band and 0.001 for H— band). 
For the brightest sources (J and H < 13.5), however, 
they show offsets of ~0.05 for both filters. Since the halo 
region near a very bright star has a dark background, 
it causes the flux of the bright star to be overestimated. 
The dark halos are induced by the sky subtraction during 
preprocessing because bright stars in the dithered images 
overestimate the sky values. The scatter of the magni¬ 
tude difference is larger for fainter sources, which may be 
caused by the photometric uncertainties for the sources. 
We calculated the variances of the offsets for both filters 
using the standard deviation of the offsets. They have 
values of ^0.05, which coincide with the zeropoint un¬ 
certainties. Therefore the magnitudes of our catalog are 
accurate to ~0.05 mag. 

4. CATALOGS PROPERTIES 
4.1. Object Detection and Photometry 

To detect sources a nd conduct photometry we em¬ 
ployed SExtractor (SE: iBertin fe Arnoutslll996l) . 

We found optimal parameters for SE to obtain high 
completeness while not sacrificing reliability of the 
detection (Sections 14.61 and 14.71) for all magnitude 
bins. To do so, we made artificial objects using the 
IRAF noao.artdata package and they are placed 
randomly at the background regions of the images. 
We performed detections and photometry with various 
parameter settings. After many trials with varying 
SE parameters, we chose the following set of detec¬ 
tion parameters for SE: DETECT_MINAREA=20, 
DETECT_THRESH=1.1, DEBLEND_NTHRESH=64, 
DEBLEND_MINCONT=0.005, BACK_SIZE=500, 
BACK_FILTERSIZE=2. The adopted detection 
parameters corresponds to 5 /N = 5. 

We detected 208,020 sources and 203,832 sources for 
J and H— band, respectively, at the magnitude limits 
brighter the 99% reliability (see Section I4~T1) . Among 
these objects, 28.3% in J and 30.8% in H are detected 
only in a single filter image. Most of the sources detected 
only in a single filter image are those at near the magni¬ 
tude limits, suggesting that they are missed in an image 
in the other band due to the detection limit. The reason 
for the percentage of sources detected only in J—band 
being sli ght ly hig her is the deep J—band catalog from 
llmai et all (120071 1. If a bright object of J < 19 mag or 
H < 20 mag has no counterpart in the other band, we 


find that they are a crosstalk or cosmic ray. We d erived 
auto magnitudes with Kron-like elliptical aperture dKronl 
I1980T) . which are assumed to be total magnitudes. For 
the aperture magnitudes, we set the aperture diameter 
size to 3 times the FWHM. We found a diameter of 3 
times the FWHM to be an optimal diameter which can 
contain the majority of the flux from an object but is not 
affected by the local fluctuation in the sky background. 
This has been checked by growth curves of objects with 
various magnitudes. The magnitude errors are computed 
by combining the standard SE errors based on Poisson 
statistics and the zeropoint errors. At the bright end, 
the photometry is free from saturation above J ~10 mag 
and H ~11 mag. 

Since the stacking of the dithered frames with sub¬ 
pixels during preprocessing causes poor measurement of 
pixel-by-pixel rms for the correlated neighboring pix¬ 
els, the errors from SE, which only consider the Pois¬ 
son noise, can be underestimated. To improve the ac¬ 
cu racy of the est i mated err ors, we adopted th e method 
in lGawiser et all (1200611 and lJeon et all (I2010T) . calculat¬ 
ing the noise rate induced by the pixel-by-pixel corre¬ 
lated noise. We place ~3000 apertures of a given size 
randomly located in the background area. The aper¬ 
ture has an area of n P i X and an effective diameter of 

N = n pix . The background rms fluctuations are derived 
from a Gaussian curve fit to the histogram of fluxes in¬ 
side each aperture. The noise calculated from the rms 
fluctuations for a given effective diameter is expressed as 
o n oc N@. If the noise is only dependent on the Poisson 
noise, it is independent from one pixel to another. There¬ 
fore for this case, on oc N. However, for the adjoining 
pixels, the sky noise has /3 > 1. For our data, we find 
on oc N 1a . For an aperture with 5.4" diameter (3 times 
the typical FWHM of our images), there are about 250 
pixels (or N ~ 16). Therefore, the photometric errors of 
faint objects (J > 18 mag) may be underestimated by a 
maximum factor of 3. For the brighter objects (J < 16 
mag), the photometric errors are not affected by this ef¬ 
fect, since their errors are dominated by noises from the 
source. 

The detection and photometry of the sources are car¬ 
ried out for each subfield image, rather than a single 
mosaic of all the subfields, since there is a variation in 
the weather condition among subfields. Catalogs from 
individual subfields are then merged together to form a 
master catalog for each filter. However, since each sub- 
field image overlaps with a neighboring images, catalogs 
created from each subfield image contain duplicate en¬ 
tries in the overlap region. To deal with the duplicate 
entries during the merging process of all the subfield cat¬ 
alogs, we set the criterion for choosing the entry with the 
smallest FLAGS values from SE and smallest photomet¬ 
ric errors for multiple entries with positional offsets less 
than 1.5". Only the object chosen with the above crite¬ 
rion, and its SE output values derived from the chosen 
subheld are included in the final catalog. 

When we merged the J and H —band catalogs, we used 
a 1.5" radius. Every row, regardless of whether they have 
information in both bands, is included in the final cata¬ 
log. Regarding the stellarity from SE, we chose a single 
stellarity value for the final catalog where the photomet¬ 
ric error of the source is smaller between the two filters. 
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4.2. Crosstalk Flagging 

The FLAMINGOS chip exhibits crosstalk artifacts 
between amplifiers on both sides of a bright source 9 , 
where the crosstalk signals can be identified as doughnut¬ 
shaped signals at every 128 pixels, which corresponds to 
the width of one amplifier, from the position of a bright 
source. The strength of the crosstalk signal is stronger 
for a brighter source, and it decreases as the distance 
from the bright source increases. We examined the mag¬ 
nitude of the source that generates crosstalks and the 
positions of the crosstalks. The crosstalks start to be 
seen on both sides of a bright source at a magnitude of 
J ~ 15.0 mag and H ~ 15.5 mag in the stacked image. 
For sources that are about 3 magnitudes brighter than 
that limit, the crosstalks continue to the end of a quad¬ 
rant. The bright sources create the crosstalk signal at 
distances of multiples of 256 pixels in our stacked image 
whose pixel sizes are a half of the pixel size of the raw 
image. The direction in which the crosstalks appeared 
depends on the readout direction and where the object is 
located, but it is predictable. Using these properties, we 
predicted the positions of the crosstalks and included the 
flag information for objects at these points in our merged 
catalog that the object is a potential crosstalk signal or 
its photometry is affected by the crosstalk. 

4.3. Data from llmai et al . I [2001 1 ) 

The J and K s —band imaging was conducted in 2004, 
covering 75 0 arcm in 2 near the center of NEP-Deep by 
llmai et all (1200711 . where we avoided duplicate observa¬ 
tion in J—band. To construct a continuous area coverage 
of the J—band imaging in th e NEP field, we included 
the J—band images from llmai et ahl (2007) when cre¬ 
ating a catalog of J—band sources. For consistency in 
the final catalog, we conducted new detections and pho¬ 
tometry using our own method (see Section 14.11) . Their 
flux cali bration was done using standard stars listed in 
iPersson et al.! (Il998l f. resulting in uncertainties in the 
flux calibration up to 0.1 mag. We were able to reduce 
the errors to less than 0.05 mag using our flux calibra¬ 
tion method (see Section m - We added this catalog 
to the final catalog to fill the gap in the NEP field from 
our J_—band survey. The average depths of images of 
llmai et al . ( 1200 711 are ~0.3 mag deeper than those of our 
data. 

4.4. Star-Galaxy Separation 

We used the same method as I Jeon et al.1 (I2010H to sep¬ 
arate galaxies from stars. Since most stars have similar 
slopes at wavelengths greater than 1/rm on their SEDs, 
NIR colors of stars lie in narrow range while those of 
distant galaxies are more diverse due to redshifted SED. 
Figure [7] shows the B — H versus H — N 2 color-color 
diagram. The stellar locus can clearly be seen by the 
clouded gray points and overlapping contour (H — N2 ~ 
0.5). Also the redshift tracks of two types of galaxies are 
shown to check the location of galaxies. The dashed line 
is a redshift track of a star-forming galaxy (M51) and 
the dotted line is that of a passively evolving 5 Gyr-old 
galaxy with a spontaneous burst, a metallicity Z=0.02, 
and the Salpeter initial mass function (IMF) created with 

9 http://flamingos.astro.ufl.edu/DetQE/index.html 


the iBruzual fc Chariot! (1200311 model. As indicated in 
Figure [71 we find that the galaxies and stars can be sepa¬ 
rated effectively with a color cut at H — N 2 = —0.2 mag, 
where objects with H — N2 < —0.2 are mostly stars. 
However, galaxies at z ~ 0 also have H — N 2 < —0.2 
mag, and there are also stars with peculiar colors with 
H — N 2 > —0.2. To deal with such case, we also use 
the stellarity parameter from SE of objects in th e op¬ 
tical catalog dJeon et alJl2010t iHwang et al.ll2007ll since 
the optical images offer better spatial resolution (~ 1"). 
The SE stellarity parameter has a value between 1 (point 
source) and 0 (extended source). For the galaxies with 
H — N 2 < —0.2 mag, they can be selected using the 
criterion stellarity < 0.8. Likewise, point sources with 
stellarity > 0.8 with H — N 2 > —0.2 mag are classi¬ 
fied as stars. This method was applied only to objects 
brighter than the 50% completeness magnitude, which 
correspond s to R < 20 mag or r < 22 (see Figure 6 in 
iJeon et all (I2010H and Figure 2 in IHwang et al.1 (I2007D L 
since the 50% completeness magnitude is the limit for 
which the stellarity parameter can be used to separate 
the extended sources from the point sources. Near the 
detection limit, all objects fainter than the magnitude of 
50% completeness are considered as galaxies. 

4.5. Catalog Format 

Table 0 shows a part of the J and H— band merged 
catalog. The catalog contains objects whose detection 
reliability is higher than 99%. The magnitude cut was 
computed independently for each field in each filter (see 
Section l4~71) . The catalog contains 15 columns; identifi¬ 
cation number (ID), the right ascension (R.A.) and dec¬ 
lination (Deck), J and if—band magnitudes and errors, 
stellarity, galactic extinction values for each filter, and 
the information about the crosstalk. The object IDs con¬ 
sist of their field name and number in the order of their 
R.A. for each field. For both J and H— band, we provide 
two kinds of magnitudes; auto magnitude (J auto and 
Hauto) using Kron-like elliptical aperture and aperture 
magnitude ( Japer and Haper) within a circular aper¬ 
ture with a diameter of three times FWHM. The mag¬ 
nitude errors are computed from the square root sum of 
the Poisson noise from the source and the zeropoint error 
from the photometric calibration. The zeropoint errors 
have values less than 0.05 mag. The distributions of the 
errors are shown in Figure[5] The large dispersions of the 
errors are induced by the different depths of each field. 
The stellarities have values between 0 (extended source) 
and 1 (point source), and the distributions of the stel¬ 
larities are also shown in Figure [8] The galactic jnrtinc- 
tion v alues are from the extinction map of Schlegcl et al. 
(119981 ). Since the color excess E(B — V) has values of 
0.034 - 0.058, Aj ~ 0.04 and Ah ~ 0.03, which are both 
smaller than the photometric zeropoint errors. Note that 
the magnitudes in the final catalog are not corrected for 
the galactic extinction. The final column contains in¬ 
formation regarding the crosstalks; the erf means that 
the object is located on a crosstalk position. We used a 
dummy value of 99.999 for non-detections. If the object 
has only single filter data, we used NaN for the other 
filter information. 

Using the filter response curve data of FLAMINGOS 10 

10 http://flamingos.astro.ufl.edu/Filter_Info/index.html 
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Figure 7. B — H versus H — N2 color-color diagram and the density contours of sources detected in B , H, and AT2—band. The points 
and the contours are made with objects that have B —band magnitude errors less than 0.1 mag. The dashed line is a reds hift track of 
a star-forming galaxy (M51) and the dotted line is for a passively evolving galaxy made with the IBruzual fe Chariot! (1200311 model of a 
passively evolving 5 Gyr-old galaxy with a spontaneous burst, a metallicity Z=0.02, and the Salpeter IMF. We use the vertical line as a 
boundary that distinguishes stellar sources against galaxies. 


and the Vega spectrum from iBohlin fe Gillilandl (l2004lh 
we computed the conversion factor between AB and Vega 
magnitude. The conversion formulae from Vega to AB 
magnitudes are as follows: 

JiVega) = J(AB) - 0.926 

H(Vega) = H{AB) - 1.367 

K s (Vega) = K S (AB) - 1.798 

We apply this method to the 2MASS filter system 
( J 2 MASS and H 2 MASS ) and found that they have val¬ 
ues similar to those of the FLAMINGOS filter system, 
0.910 for J 2 mass and 1.383 for H 2 mass- 

4.6. Completeness 

To estimate the completeness for each field, we used ar¬ 
tificial objects from the noao. artdata package of IRAF. 


The artificial objects are created with magnitude be¬ 
tween 11 and 25 and they are placed randomly in the 
background regions of the images. The number of 
point sources and extended sources created are 100 each. 
Among the extended sources, 40% are elliptical galax¬ 
ies while the others are disk galaxies, with a minimum 
redshift of z ~ 0.01. The angular sizes of the galaxies, 
which are equivalent to the half-flux diameters of the 
model objects, are set to be the FWHM values for each 
image. We performed detections and photometry with 
the same parameter settings of SE which were applied 
for the actual images. We define the completeness as 
the ratio of the number of detected artificial sources to 
the number of input objects in a given magnitude bin. 
Figure [9] represents the completeness of J and H —band 
data as a function of magnitude. The gray lines repre¬ 
sent the completeness for each subfield. The thick, black 
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The photometric error distribution and the stellarity distribution 


Figure 8. 

lines are for the deepest (FLM27 in J and FLM58 in 
H) and the shallowest (FLM19 in J and FLM07 in H) 
fields according to the 5er detection limits. In Figure 0 ] 
we compared the magnitude of 50% completeness (aster¬ 
isks) for each field with the 5er detection limit, the 99% 
reliability (see Section 14771 ) . the Z p , and the total expo¬ 
sure time. The magnitudes of 50% completeness for each 
field are listed in Table El We can see that the magni¬ 
tudes of 50% completeness are ~1 mag brighter than the 
5cr detection limits. 

4.7. Reliability 

To examine the reliability of our source detection, we 
used the false detection rate for each field and filter. The 
false detection rate is defined as the ratio of the number of 
objects detected in a negative image to the number of ob¬ 
jects detected in the original image in a given magnitude 
bin. The negative images were created by multiplying 
—1 to the original images. If an image has only noise, its 
false detection rate would be 1. We performed detections 
and photometry on the negative image with the same pa¬ 


rameter settings of SE which were applied for the actual 
images. Figure [TU] represents the false detection rates 
of J and IF—band data as a function of magnitude. The 
gray lines represent the false detection rates for each sub¬ 
field. The thick, black lines are for the deepest (FLM27 
in J and FLM58 in H) and the shallowest (FLM19 in J 
and FLM07 in H) fields according to the 5cr detection 
limits. In Figure 01 we compared the magnitude of 1% 
false detection rates (equivalent to 99% reliability) (di¬ 
amonds) for each field with the 5 a detection limit, 50% 
completeness, the Z pi and the total exposure time. The 
magnitude of 99% reliability for each field are presented 
in Table [G] We can see that the magnitudes of 99% reli¬ 
ability are similar with the 5cr detection limits. 

5. PROPERTIES OF THE DATA 
5.1. Source Number Counts 

The galaxy number count as a function of magnitude is 
a quantity used in investigating galaxy evolution. Due to 
the contamination of dust extinction and star formation, 
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Figure 9. Completeness (the ratio of the number of detected artificial sources to the number of input objects in a given magnitude bin) 
for J and H— band. The gray solid lines are the completeness of all subfields from the simulation method, especially for the deepest fields 
(FLM27 in J and FLM58 in H; the upper black solid lines) and the shallowest fields (FLM19 in J and FLM07 in H\ the lower black solid 
lines) in the 5<r detection limits. The horizontal dashed line is for 50% completeness cut. The magnitudes of 50% completeness for each 
subfield are shown in Table [6] 
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Figure 10. False detection rate (the ratio of the number of objects detected in the negative image to the number of objects detected in 
the original image for each magnitude bin) of all subfields in J, and H— band. The black solid lines are for the deepest or the shallowest 
fields as same as Figure [9] The magnitudes of false detection rate 1% for each field are shown in Table [6] 




































NEP FLAMINGOS NIR Catalog 


13 




16 17 18 19 20 21 22 

H magnitude [AB] 


Figure 11. Galaxy (two solid lines with crosses) and stellar (solid lines with asterisks) number counts of J and H— band. The error 
estimates of number counts are from the Poisson errors. The lower solid line is for galaxy number count without completeness correction, 
and the upper solid line is from the galaxy number count corrected using the completeness. We compare our number counts with the 
number counts in the literature (other symbols). The galaxy and stellar number count data are provided in Table |4] 
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Table 4 

Galaxy and stellar number counts. 


mag 


J- 

band 



H- 

-band 


Galaxy 

Number 

Error 

Corrected 

Number 

Star 

Number 

Galaxy 

Number 

Error 

Corrected 

Number 

Star 

Number 

16.25 

4.05e+01 

1.98e+01 

9.94e-01 

8.31e+02 

4.45e+01 

2.05e+01 

9.94e-01 

8.31e+02 

16.75 

6.59e+01 

2.54e+01 

9.85e-01 

1.02e+03 

7.89e+01 

2.78e+01 

9.87e-01 

1.02e+03 

17.25 

1.02e+02 

3.15e+01 

9.78e-01 

1.27e+03 

1.35e+02 

3.64e+01 

9.78e-01 

1.27e+03 

17.75 

1.82e+02 

4.23e+01 

9.67e-01 

1.54e+03 

2.47e+02 

4.98e+01 

9.63e-01 

1.54e+03 

18.25 

3.13e+02 

5.59e+01 

9.52e-01 

1.84e+03 

4.04e+02 

6.37e+01 

9.37e-01 

1.84e+03 

18.75 

4.67e+02 

6.84e+01 

9.25e-01 

2.18e+03 

6.28e+02 

7.99e+01 

8.90e-01 

2.18e+03 

19.25 

7.06e+02 

8.45e+01 

8.71e-01 

2.97e+03 

8.96e+02 

9.54e+01 

7.94e-01 

2.97e+03 

19.75 

1.01e+03 

1.01e+02 

7.73e-01 

4.27e+03 

1.31e+03 

1.15e+02 

6.42e-01 

4.27e+03 

20.25 

1.64e+03 

1.27e+02 

6.17e-01 

5.89e+03 

3.00e+03 

1.70e+02 

4.42e-01 

5.89e+03 

20.75 

4.96e+03 

2.23e+02 

4.54e-01 

7.50e+03 

8.11e+03 

2.87e+02 

2.52e-01 

7.50e+03 


Note. — Units of the galaxy and stellar number densities and their errors are numbers 0.5 mag 
deg” 2 . 


the galaxy number counts in NIR bands provide more 
accurate information of galaxy evolution than those using 
the optical bands. Here, we show the galaxy number 
count as a function of magnitude for representation of 
depth and homogeneity of our survey. The total area 
used in calculating the ga laxy number count is 4.98 deg 2 , 
where the optical data dJeon et al.1 120101 IHwang et alj 
l2007f ) and the AKA HI III data overlap. 

Using the galaxies defined so in Section 14.41 we cal¬ 
culated the number of galaxies per square degree for 
each magnitude bin. Since the depth of each subfield is 
different from another, we computed the number count 
for each subfield down to the 99% reliability magnitudes 
limit. Then the completeness correction was applied to 
the number counts for each subfield. Finally, we com¬ 
bined the number counts from all the subfields to derive 
the final number count result. We also calculated the 
number count of stars up to the 99% reliability magni¬ 
tude. The errors from the number counts are computed 
from the Poisson errors. 

Figure [IT] shows the galaxy and star number counts 
in J and H— band. The lower solid line with crosses 
represents the raw number counts and the upper 
solid line with crosses shows the number counts af¬ 
ter applying the completeness correction. We com¬ 
pared our number counts with those in the literature 


(lYan et al. 

1998): 

Saracco et al.l 

1999, 

2QQl| Martini 

2001 

Chen et al. 

2002 

: Iovino et al. 

2005 

; Frith et al.| 

2006; 

Metcalfe et al. 2009: Cristobal-Hornillos et al.l 

2009; 

Keenan et al. 20101. Since the filter systems in the liter- 


ature do not exactly match with that of our survey and 
also due to th e cosmic varian ce, so me offsets may exist. 
For example, iMetcalfe et al.l (120061 ) mention that their 
number counts, covering only 0.95x0.95 arcmin 2 , can be 
uncertain by up to a factor of 2 due to the cosmic vari¬ 
ance. We can see that our number counts are consistent 
with other results. The solid line with asterisks shows 
the number counts of stars. The galaxy and star number 
count data are provided in Table [4] 

5.2. Color-Magnitude Diagram 

Figure [T21 shows the color-magnitude diagram of J ver¬ 
sus J — H and its histogram. The upper panel shows the 
distribution for galaxies with J—band magnitude error 
less than 0.1 mag, and the lower panel shows that for 


stars with the same condition as galaxies in the upper 
panel. We used the redshift track to check the posi¬ 
tion of galaxies on the color-magnitude diagram. We 
used the same conditions for the galaxy redshift track 
as used in Section 14.41 These redshift tracks represent 
galaxies with characteristic absolute magnitude M*. We 
adopt M* — 51og fo for J—band as —22.85 ± 0.04 from 
I Jones et~ahl (l2006i l. We c alculated the fc-co rrection val¬ 
ues in a standard manner ([Hogg et al .1120021) . The J — H 
colors of the star forming galaxies and passive evolv¬ 
ing galaxies are similar to each other up to redshift ~1. 
Therefore the galaxies on Figure [12] are those within the 
redshift range 0 to 1, and those at higher redshift with 
magnitude brighter than M*. 

6. SUMMARY 

We have presented the properties of sources using J 
and H —band data from FLAMINGOS on the KPNO 
2.1m telescope, a set of follow-up imaging observations of 
the AKARI NEP field. We cover 5.1 deg 2 with 5cr depth 
of ~21.6 mag and ~21.3 mag (AB) for J and H— band, 
respectively. We detected 208,020 sources for J—band 
and 203,830 sources for H— band. Our NIR data provides 
multi-wavelength coverage for the AKARI NEP sources 
along with the optical and the MIR data, enabling de¬ 
tailed SED analysis. Through these SEDs, we can obtain 
key physical properties such as stellar mass and photo¬ 
metric redshifts, and select interesting objects such as 
high redshift galaxies. 
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Table 5 

Observation summary. 


Filter 

Total exposure time [minutes] 

Exposure time per image [minutes] 

Seeing ["] 

5 (j detection limits [AB mag] 

50% Completeness [AB mag] 

J 

28.5 - 52.5 (36) 

1.5 

1.4 - 2.4 (1.8) 

20.88 - 21.94 (21.63) 

20.11 21.01 (20.65) 

H 

11.3 - 34.0 (21.3) 

0.5 - 1.0 (0.7) 

1.4 - 2.6 (1.7) 

20.42 - 21.72 (21.29) 

19.44 - 20.81 (20.35) 


Note. 


The values in the parentheses are the median values. 



Table 6 

Photometric zeropoint, seeing, 5<r detection limit, magnitude of 50% 
completeness, and magnitude of 99% reliability for each subfield. 


Field 



J —band 





H —band 



Z P a 

(mag) 

Seeing 

(") 

detection limits b 
(mag) 

comp50% c 

(mag) 

relia99% d 

(mag) 

(mag) 

Seeing 

(") 

detection limits b 
(mag) 

comp50%° 

(mag) 

relia99% d 

(mag) 

FLM01 

22.43 

1.67 

21.31 

20.40 

21.46 

23.22 

1.98 

21.06 

20.06 

21.12 

FLM02 

22.41 

1.67 

21.87 

20.85 

21.78 

23.26 

2.43 

20.93 

19.85 

21.15 

FLM03 

22.43 

1.67 

21.85 

20.97 

21.74 

23.26 

2.58 

20.74 

19.70 

21.11 

FLM04 

22.42 

1.67 

21.70 

20.69 

21.63 

23.22 

1.67 

21.17 

20.18 

21.34 

FLM05 

22.45 

1.98 

21.24 

20.20 

21.30 

23.21 

1.67 

21.08 

20.14 

21.12 

FLM06 

22.43 

1.98 

21.32 

20.25 

21.47 

23.14 

2.28 

20.49 

19.42 

21.01 

FLM07 

22.36 

1.98 

21.49 

20.47 

21.48 

23.16 

1.98 

20.42 

19.58 

20.63 

FLM08 

22.46 

1.98 

21.48 

20.49 

21.48 

23.13 

1.67 

20.96 

19.94 

20.94 

FLM09 

22.41 

1.67 

21.74 

20.78 

21.58 

23.34 

1.82 

20.92 

19.97 

20.94 

FLM10 

22.43 

1.67 

21.73 

20.80 

21.68 

23.21 

1.98 

20.94 

20.08 

21.11 

FLM11 

22.42 

1.98 

21.17 

20.09 

21.33 

23.23 

1.52 

21.28 

20.35 

21.20 

FLM12 

22.57 

1.67 

21.72 

20.79 

21.54 

23.28 

1.52 

21.24 

20.24 

21.14 

FLM13 

22.41 

1.82 

21.71 

20.69 

21.69 

23.22 

1.52 

21.30 

20.37 

21.24 

FLM14 

22.71 

1.82 

21.61 

20.70 

21.65 

23.22 

1.52 

21.31 

20.40 

21.29 

FLM15 

22.40 

1.82 

21.76 

20.70 

21.75 

23.21 

1.67 

21.29 

20.29 

21.16 

FLM16 

22.58 

1.98 

21.53 

20.38 

21.64 

23.27 

1.67 

21.00 

20.02 

21.06 

FLM17 

22.48 

1.82 

21.58 

20.66 

21.61 

23.19 

1.82 

21.01 

19.98 

21.09 

FLM18 

22.49 

1.98 

21.34 

20.31 

21.43 

23.26 

1.98 

20.91 

19.97 

21.08 

FLM19 

22.44 

2.13 

21.06 

20.10 

21.25 

23.25 

1.67 

21.28 

20.33 

21.17 

FLM20 

21.98 

1.82 

21.47 

20.36 

21.45 

23.21 

2.13 

21.17 

20.10 

21.33 

FLM21 

22.47 

1.52 

21.63 

20.73 

21.58 

23.22 

1.98 

21.28 

20.20 

21.14 

FLM22 

22.51 

1.52 

21.63 

20.62 

21.51 

23.27 

1.52 

21.47 

20.46 

21.28 

FLM23 

22.47 

1.82 

21.93 

20.98 

21.69 

23.23 

1.52 

21.48 

20.58 

21.21 

FLM24 

22.52 

2.13 

21.68 

20.59 

21.76 

23.24 

1.67 

21.36 

20.31 

21.27 

FLM26 

22.59 

2.43 

21.51 

20.26 

21.65 

22.81 

1.52 

21.70 

20.71 

21.38 

FLM27 

22.53 

1.67 

21.94 

20.93 

21.67 

23.20 

1.82 

21.68 

20.71 

21.45 

FLM28 

22.54 

1.98 

21.25 

20.22 

21.31 

23.46 

1.52 

21.30 

20.48 

21.15 

FLM29 

22.35 

1.67 

21.69 

20.70 

21.70 

23.36 

1.82 

21.16 

20.19 

21.15 

FLM30 

21.78 

1.67 

21.53 

20.34 

21.56 

23.29 

2.13 

21.12 

20.17 

21.35 

FLM31 

22.39 

1.82 

21.50 

20.58 

21.59 

23.37 

1.52 

21.38 

20.55 

21.20 

FLM32 

22.40 

1.67 

21.75 

20.73 

21.59 

23.23 

1.67 

21.29 

20.42 

21.21 

FLM33 

22.43 

1.98 

21.72 

20.78 

21.64 

23.25 

1.67 

21.45 

20.58 

21.36 

FLM34 

22.43 

1.67 

21.73 

20.70 

21.50 

22.81 

1.67 

21.20 

20.24 

21.14 

FLM35 

22.30 

2.13 

21.80 

20.50 

21.78 

22.84 

1.52 

21.67 

20.64 

21.39 

FLM36 

22.61 

2.13 

21.60 

20.53 

21.63 

23.08 

1.67 

21.23 

20.27 

21.13 

FLM37 

22.50 

1.52 

21.62 

20.67 

21.50 

23.29 

1.67 

21.22 

20.25 

21.21 

FLM38 

22.40 

1.52 

21.87 

20.80 

21.44 

23.29 

1.67 

21.25 

20.18 

21.12 

FLM39 

22.39 

1.67 

21.76 

20.72 

21.65 

23.22 

1.98 

21.19 

20.17 

21.27 

FLM40 

22.36 

1.98 

21.41 

20.50 

21.56 

23.28 

1.67 

21.29 

20.34 

21.20 

FLM41 

22.23 

1.82 

21.21 

20.14 

21.29 

23.23 

1.52 

21.33 

20.48 

21.21 

FLM42 

22.44 

1.98 

21.65 

20.66 

21.67 

23.24 

1.67 

21.32 

20.46 

21.20 

FLM43 

22.54 

1.98 

21.65 

20.61 

21.69 

22.82 

1.67 

21.28 

20.40 

21.22 

FLM44 

22.53 

1.82 

21.76 

20.71 

21.79 

22.81 

1.82 

21.14 

20.22 

21.22 

FLM45 

22.40 

1.52 

21.69 

20.77 

21.61 

22.80 

1.67 

21.48 

20.51 

21.36 

FLM46 

22.58 

2.28 

21.50 

20.22 

21.63 

22.87 

1.67 

21.41 

20.42 

21.22 

FLM47 

22.38 

1.82 

21.38 

20.41 

21.45 

23.27 

1.67 

21.33 

20.34 

21.23 

FLM48 

22.41 

1.82 

21.41 

20.37 

21.52 

23.21 

1.67 

21.35 

20.34 

21.38 

FLM49 

22.41 

1.67 

21.55 

20.47 

21.44 

23.27 

1.67 

21.37 

20.32 

21.34 

FLM50 

22.59 

1.67 

21.47 

20.51 

21.50 

23.20 

1.67 

21.42 

20.42 

21.40 

FLM51 

22.42 

1.67 

21.66 

20.76 

21.65 

23.19 

1.82 

21.38 

20.40 

21.40 

FLM52 

22.41 

1.67 

21.72 

20.60 

21.59 

23.32 

1.67 

21.34 

20.23 

21.25 

FLM53 

22.56 

1.67 

21.70 

20.86 

21.68 

23.17 

1.37 

21.64 

20.80 

21.27 
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Table 6 


Continued 


Field 



J —band 





H —band 



(mag) 

Seeing 

(") 

detection limits b 
(mag) 

comp50% c 

(mag) 

relia99% d 

(mag) 

(mag) 

Seeing 

(") 

detection limits b 
(mag) 

comp50%° 

(mag) 

relia99% d 

(mag) 

FLM54 

22.56 

1.67 

21.65 

20.62 

21.57 

23.24 

1.52 

21.42 

20.54 

21.22 

FLM55 

22.46 

1.67 

21.80 

20.75 

21.66 

23.28 

1.67 

21.37 

20.40 

21.30 

FLM56 

22.51 

1.82 

21.68 

20.45 

21.79 

23.19 

1.67 

21.28 

20.25 

21.40 

FLM57 

22.49 

1.98 

21.50 

20.44 

21.52 

23.19 

1.67 

21.36 

20.47 

21.22 

FLM58 

22.50 

1.67 

21.51 

20.43 

21.50 

23.20 

1.37 

21.72 

20.75 

21.38 

FLM59 

22.48 

1.67 

21.59 

20.63 

21.53 

23.20 

1.52 

21.57 

20.61 

21.40 

FLM60 

22.49 

1.67 

21.67 

20.68 

21.65 

23.20 

1.52 

21.59 

20.59 

21.36 

FLM61 

22.56 

1.67 

21.75 

20.68 

21.57 

23.86 

1.52 

21.03 

20.19 

20.96 

FLM62 

22.62 

1.67 

21.81 

20.74 

21.71 

23.32 

1.52 

21.48 

20.60 

21.21 

FLM63 

22.52 

1.67 

21.81 

20.78 

21.48 

23.18 

1.52 

21.50 

20.64 

21.26 

FLM64 






23.16 

1.82 

21.25 

20.31 

21.20 

FLM65 






22.79 

1.82 

21.10 

20.12 

21.10 

FLM66 






22.83 

1.67 

21.32 

20.49 

21.23 


a The apparent magnitude is M = Z p — 2.5 log ( DN/sec ). 

k The 5 a detection limit for a circular aperture with diameter of 3 times FWHM 
c The magnitude corresponding to the 50% completeness 
^ The magnitude corresponding to the 99% reliability 
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Table 7 

J and H —band merged catalog. 


to 

o 


ID 

(1) 

R.A. 

(2) 

Decl. 

(3) 

Jauto 

(4) 

Error 

(5) 

Japer 

(6) 

Error 

(7) 

a AUTO 
(8) 

Error 

(9) 

Haper 

(10) 

Error 

(ii) 

Stellarity 

(12) 

Aj 

(13) 

A h 

(14) 

flag 

(15) 

FLM01.00651 

18:03:16.74 

+67:32:13.20 

20.911 

0.136 

20.966 

0.168 

20.885 

0.177 

21.299 

0.276 

0.65 

0.045 

0.029 


FLM01.00652 

18:03:16.75 

+67:41:18.77 

NaN 

NaN 

NaN 

NaN 

20.838 

0.202 

20.746 

0.172 

0.47 

NaN 

0.027 


FLM01.00653 

18:03:16.81 

+67:22:15.77 

21.172 

0.199 

21.199 

0.200 

99.999 

99.999 

99.999 

99.999 

0.48 

0.045 

NaN 


FLM01.00654 

18:03:16.88 

+67:37:10.82 

20.198 

0.074 

20.109 

0.092 

99.999 

99.999 

99.999 

99.999 

0.67 

0.042 

NaN 


FLM01.00655 

18:03:16.91 

+67:36:12.73 

99.999 

99.999 

99.999 

99.999 

20.715 

0.158 

20.921 

0.199 

0.49 

NaN 

0.027 


FLM01.00656 

18:03:16.92 

+67:30:57.40 

20.254 

0.086 

20.342 

0.106 

20.400 

0.112 

20.464 

0.137 

0.89 

0.046 

0.029 


FLM01.00657 

18:03:16.96 

+67:33:30.66 

20.399 

0.103 

20.524 

0.120 

19.795 

0.100 

19.798 

0.087 

0.16 

0.044 

0.028 


FLM01.00658 

18:03:16.97 

+67:26:35.22 

20.136 

0.072 

20.265 

0.101 

99.999 

99.999 

99.999 

99.999 

0.90 

0.046 

NaN 


FLM01.00659 

18:03:17.05 

+67:27:24.29 

99.999 

99.999 

99.999 

99.999 

21.049 

0.177 

21.198 

0.253 

0.61 

NaN 

0.030 


FLM01.00660 

18:03:17.05 

+67:37:27.81 

20.610 

0.184 

20.866 

0.155 

99.999 

99.999 

99.999 

99.999 

0.55 

0.042 

NaN 


FLM01.00661 

18:03:17.09 

+67:36:35.88 

21.122 

0.187 

21.149 

0.195 

99.999 

99.999 

99.999 

99.999 

0.58 

0.042 

NaN 


FLM01.00662 

18:03:17.10 

+67:37:22.53 

21.117 

0.217 

21.255 

0.214 

99.999 

99.999 

99.999 

99.999 

0.65 

0.042 

NaN 


FLM01.00663 

18:03:17.15 

+67:38:35.97 

20.786 

0.155 

20.829 

0.151 

99.999 

99.999 

99.999 

99.999 

0.51 

0.041 

NaN 


FLM01.00664 

18:03:17.17 

+67:29:08.37 

19.167 

0.070 

19.209 

0.066 

19.016 

0.067 

19.050 

0.065 

0.89 

0.046 

0.030 


FLM01.00665 

18:03:17.19 

+67:41:50.41 

NaN 

NaN 

NaN 

NaN 

18.020 

0.058 

18.138 

0.057 

0.17 

NaN 

0.027 


FLM01.00666 

18:03:17.22 

+67:34:37.53 

20.523 

0.148 

20.908 

0.160 

99.999 

99.999 

99.999 

99.999 

0.34 

0.043 

NaN 


FLM01.00667 

18:03:17.23 

+67:32:17.26 

20.957 

0.135 

21.282 

0.219 

20.929 

0.197 

20.760 

0.174 

0.40 

0.045 

0.029 


FLM01.00668 

18:03:17.23 

+67:33:04.26 

20.589 

0.161 

20.984 

0.170 

99.999 

99.999 

99.999 

99.999 

0.47 

0.044 

NaN 


FLM01.00669 

18:03:17.24 

+67:37:29.69 

99.999 

99.999 

99.999 

99.999 

19.754 

0.106 

20.266 

0.118 

0.49 

NaN 

0.026 


FLM01.00670 

18:03:17.26 

+67:39:22.08 

21.222 

0.188 

21.120 

0.191 

99.999 

99.999 

99.999 

99.999 

0.65 

0.042 

NaN 


FLM01.00671 

18:03:17.31 

+67:45:01.33 

NaN 

NaN 

NaN 

NaN 

20.456 

0.114 

20.455 

0.136 

0.37 

NaN 

0.028 


FLM01.00672 

18:03:17.35 

+67:40:17.33 

21.132 

0.146 

21.319 

0.226 

99.999 

99.999 

99.999 

99.999 

0.49 

0.042 

NaN 


FLM01.00673 

18:03:17.42 

+67:35:23.47 

21.448 

0.142 

23.891 

2.331 

99.999 

99.999 

99.999 

99.999 

0.52 

0.042 

NaN 


FLM01.00674 

18:03:17.52 

+67:23:31.96 

20.646 

0.077 

20.710 

0.137 

99.999 

99.999 

99.999 

99.999 

0.70 

0.045 

NaN 


FLM01.00675 

18:03:17.58 

+67:34:33.18 

20.243 

0.142 

20.824 

0.150 

99.999 

99.999 

99.999 

99.999 

0.65 

0.043 

NaN 


FLM01.00676 

18:03:17.60 

+67:35:51.45 

19.728 

0.079 

20.022 

0.088 

99.999 

99.999 

99.999 

99.999 

0.97 

0.042 

NaN 

crt 

FLM01.00677 

18:03:17.61 

+67:30:59.62 

99.999 

99.999 

99.999 

99.999 

21.111 

0.144 

22.253 

0.654 

0.49 

NaN 

0.029 


FLM01.00678 

18:03:17.62 

+67:33:00.44 

99.999 

99.999 

99.999 

99.999 

21.034 

0.162 

21.097 

0.231 

0.54 

NaN 

0.028 


FLM01.00679 

18:03:17.63 

+67:30:23.69 

21.848 

0.234 

22.219 

0.503 

20.091 

0.111 

20.073 

0.103 

0.52 

0.046 

0.029 


FLM01.00680 

18:03:17.67 

+67:39:53.53 

17.777 

0.061 

18.106 

0.059 

17.814 

0.058 

17.948 

0.056 

0.30 

0.042 

0.027 



Note. — (1) The identification number of the source (ID). Each ID consists of the subfield name (Figure [2j and number in the order of their R.A.. 

(2) and (3) The J2000.0 R.A. and Decl. of a source in a sexigesimal. 

(4)—(11) are the total magnitude, its magnitude error, aperture magnitude and its magnitude error of a source in the J and H— band. The aperture magnitudes are obtained with 
the aperture diameter of 3 times FWHM. The magnitudes are not corrected for the galactic extinctions. We used a dummy value of 99.999 for non-detections. If the object has only 
single filter data, we used NaN for the other filter information. 

(12) The stellarity index of a source from SE. See Figure [8] for the distribution o f this i n dex. _ 

(13-14) The galactic extinction values of J and H —band from extinction map o fISchlegel et ahl d 19981) . 

(15) The flag parameter. It contains the crosstalk information ( crt ). See section T4. 21 for details . 
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